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Abstract 

A trial application of the method of Feynman-Kleinert approximants is made 
to perturbation series arising in connection with the lattice Schwinger model. 
In extrapolating the lattice strong-coupling series to the weak-coupling contin- 
uum limit, the approximants do not converge well. In interpolating between 
the continuum perturbation series at large fermion mass and small fermion 
mass, however, the approximants do give good results. In the course of the 
calculations, we picked up and rectified an error in an earlier derivation of the 
continuum series coefficients. 

PACS Indices: 11. 15. Me, 02.30.Mv, 12.20.Ds 
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I. INTRODUCTION 



The analytic continuation of a function, whose perturbation series expansion about a 
given point is known, is an important and long-standing problem. Some standard methods 
are the use of Pade approximants, or integrated differential approximants, as reviewed for 
instance by Guttmann [1]. In some circumstances, these techniques can work extremely 
well; but in others they fail. A well known difficult case is that of a function whose series 
expansion about the origin in a variable x, say, is known up to a given order, and we want to 
extrapolate to x — > oo. If there is a cut or an essential singularity at infinity, the standard 
methods do not work well. This is a common situation in lattice gauge theory, for instance. 

In recent years, a new method has been developed, involving "Feynman-Kleinert approxi- 
mants" [2]. This is a variational technique, arising originally from a Gaussian approximation 
to the path integral [3] . It has been used with great success on a number of problems, such 
as extrapolation of the weak-coupling series for the anharmonic oscillator [4], where the 
asymptotic strong-coupling behaviour was obtained with extraordinary accuracy, up to 20 
significant digits. For a discussion and further references, see the recent paper by Kleinert 
[5]. 

Here we attempt to apply the new method to a lattice gauge theory, namely the Schwinger 
model in (1+1) dimensions. In the series approach to lattice gauge theory, one generates 
strong-coupling series expansions for the physical quantities of interest, and then tries to 
extrapolate to the continuum, weak-coupling limit. It would be an important breakthrough 
if a better technique could be found for performing this extrapolation. In the report we set 
out to test whether Feynman-Kleinert approximants can fulfill this need. Unfortunately, we 
find that the technique is not very successful in this context. 

Another common problem is that of interpolation of a function, where a number of 
series coefficients are known of both the weak-coupling and the strong-coupling expansions. 
An example of this arises for the low-lying energy eigenvalues in the continuum Schwinger 
model, as functions of m/g, where m is the fermion mass and g is the electric coupling. We 
show that Feynman-Kleinert approximants can provide an effective interpolation between 
these two limits. 



II. THEORY 

For reference, let us briefly recapitulate the algorithm set out by Kleinert [2] using his 
notation. Consider a quantity E whose perturbation series expansion in some parameter a 
is known to order N: 

N 

E N = Y j a n a n . (1) 

?1=0 

Write 

N / a \ n 

E% = ^Y.^{-) , (2) 

n=0 VCJ J 

where uj is an auxiliary parameter whose value will eventually be set equal to 1, and p, q are 
indices which will be chosen later to suit the problem at hand. Now set uj = ^/Q 2 + uj 2 — Q 2 , 
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and re-expand E N in powers of A = 1 — u 2 /Q 2 , treating A as a quantity of order a, and 
truncating the re-expanded series at order N. The result is a function 



N 



W N (a,Q)=Wj2 a nfn{ty 



n=0 



a 



(3) 



where 



(4) 



Forming the first and second derivatives of Wn(cx, Q) with respect to (and setting u — 1), 
we find the positions of the turning points. The smallest among these is denoted Qn- The 
resulting Wn(cx) = Wn(cx,Qn) gives the desired approximation to the function E at finite 
a, an extrapolation of the function En to finite a. 

Now consider the limit a — > oo. At large a, fijv will scale with a like 



so that 

W N (a,Q N ) ~ 

The full 'strong-coupling' (large a) expression can now be obtained by writing 

w N (a,n) = n p w N (&,uj 2 /n 2 ), 



(5) 



(6) 



with a = a/Q q , and expanding in powers of u 2 /fl 2 , which for a 
{l/c 2 ){a/u q )- 2 l q . The result is 



W n (a) 



a 



p/q 



6d(c) + 6i(c)(-J 



a\- 2 / q 



+ . . . 



(7) 

oo behaves like 
(8) 



with 



N N-l 
/=0 j=n 



b n (c) = ±w%\a,0) « M/ % =1/c9 
(P"f?)/2^ f J 



(9) 



The indices p and q in the expansion (2) are now chosen so as to give the correct asymptotic 
behaviour (i.e. the correct powers of a) in the strong-coupling expansion (8), which we 
assume is known a priori. The leading coefficient c in (5) is found by searching for the 
extrema of the leading coefficient 6 (c) as a function of c and choosing the smallest of them. 

Finally, we have to account for the fact that Qn will have corrections to the leading 
behaviour ct^^c. We can allow for this by replacing c by 



c(a) = c + ci ( —J + c 2 



a \ 



-2/9 



a 



s -4/9 



+ 



(10) 



3 



requiring a re-expansion of the coefficients b n (c) in (9). The expansion coefficients c n are 
determined by looking for the turning points of &2n( c ), successively. 
The equations 

dhn{c) = 0, (f f) 



dc. 



n 



together with knowledge of the coefficients {a n , n = 0, . . . , N}, can in principle be used to 
determine a set of parameters {c n ,n = 0, . . . ,N}, and hence estimate the leading strong- 
coupling coefficients b n (c). This provides estimates of the asymptotic behaviour of the func- 
tion E. In the case of the strong-coupling anharmonic oscillator, it has been shown [4] 
that the resulting estimates converge exponentially fast, and have been used to estimate the 
energy eigenvalue in the strong-coupling limit accurate to 20 decimal places! 

Alternatively, any knowledge of the strong-coupling coefficients bo, b\, . . . can be fed back 
into the set of equations (9) and (11) to obtain estimates of further weak-coupling coefficients 
ajv+i, cln+2, ■ ■ ■, and thus carry the extrapolation Wn(cx) to higher orders. This is the basis 
of the interpolation algorithm between weak and strong-coupling for the function E. 

III. RESULTS 

A. Extrapolation to the Continuum Limit 

We have applied the Feynman-Kleinert approximation method to the lattice Schwinger 
model, as formulated by Banks et al. [6] . The rescaled lattice Hamiltonian is 

W = —zH = W + xV, (12) 
ag 2 



where 



Wo ^E^W + ^B- 1 )*^), (is) 

n n 

V = -iJ2 [^ t (n)e w(n V(n + 1) - h.c] , (14) 



and 



/^=^T, x= ~^l- ( 15 ) 
g z a g z a z 

Here m is the electron mass, g is the electric coupling constant, and a is the lattice spacing. 
The field <p{n) is the lattice fermion field defined on sites n, 

{0(n),0 f (m)} = 5 nm , (16) 

while L(n) is the electric field on link n, which takes only integer eigenvalues. 

Expansions of the energy eigenvalues of the operator W have been calculated to high 
orders in the perturbation parameter x, using efficient linked cluster methods [7]. Our object 
now is to see whether these series can be accurately extrapolated to the continuum limit, 
which corresponds to a — > or x — > oo, using the Feynman-Kleinert approach. 
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1. Ground State Energy 



For the ground state energy per site lvq/M, where M is the number of states, the series 
expansion is of the form 

u /M = J2a n (x 2 ) n , (17) 

n 

where a = 0; and the asymptotic behaviour as x — > oo is known [6] 

2x 

uj /M ~ . (18) 

7T 

The coefficients a n have been calculated up to order n — 15 [7]. In this case, we can apply 
the formalism of section II, taking E = ujq/M, a = x 2 , p = 2, q = 4, which implies the 
asymptotic behaviour as x — > oo 

cu /M ~ 6 x + 6i + 0(x~ 1 ). (19) 

Fig. 1 shows the resulting approximants to the ground-state energy at m/g = compared 
with numerical results obtained by other methods, and the known continuum limit, Eq. (18). 
It can be seen that the approximants do seem to converge to the correct limit, but only rather 
slowly. 

2. Energy Gap 

The energy gap to the lowest-lying positronium state has an expansion 

uji - u = a n (a; 2 ) (20) 

n 

and the expected asymptotic behaviour as x — > oo is 

ui - ujq ~ b x 1/2 . (21) 

The coefficients a n for this series have been calculated up to order n = 14 [7]. Here, we take 
a = x 2 , p = 2, q = 8, assuming the asymptotic behaviour as x — > oo 

wi - uj ~ b x 1/2 + h + 0(x~ 1/2 ) (22) 

which accords with previous studies [7,8]. In this case, however, the Feynman-Kleinert 
approximants do not converge at all at large x. Fig. 2 shows the dependence of the estimated 
value of bo as a function of the parameter Co, for the case m/g = 0, at various orders N. It can 
be seen that as Co increases, the estimate of bo comes in from infinity, oscillates around the 
correct value (marked by a solid line ), and then slowly drifts away again. Unfortunately, the 
oscillations increase with order N, and the values at the lowest turning point correspondingly 
diverge away from the correct value as iV increases. The same is true of the values at the 
lowest point of inflection. It seems that the Feynman-Kleinert approximant method fails in 
this case. 
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B. Interpolation of Continuum Series 

Another problem connected with series expansions for the Schwinger model arises in 
connection with the continuum field theory, rather than the lattice model. A few coefficients 
are known of the series expansions for the lowest-lying excited states of the model, the 
"vector" and "scalar" positronium states, in both the strong-coupling and weak-coupling 
regimes. At strong coupling, m/g — > 0, the first three coefficients have been computed by 
Adam [9], using the bosonic form of the theory: 

E (M^A^^ (23) 

9 \ 9 J n 

At weak coupling, m/g — > oo, the first three coefficients have been calculated by Sriganesh 
et al. [11], following a treatment of Coleman [10]: 

9 Vm/ V 

The problem now is, how can one compute a sensible interpolation between these two 
series to obtain accurate estimates of the energy eigenvalues at finite m/g! One possibility 
is to use "2-point" Pade approximants: but they involve only integer powers of the variable, 
{m/g) in this case. By fitting an appropriate power of the function E/g one can mimic the 
leading behaviour as m/g — > oo, but cannot reproduce the sub-leading powers. Furthermore, 
experience shows that 2-point Pade approximants are prone to exhibit spurious 'wiggles' at 
intermediate couplings between the two limits. Here we set out to apply Feynman-Kleinert 
approximants to the problem. 

Setting a = m/g, p = 1, q = 3, we have computed Feynman-Kleinert approximants 
for the energy eigenvalues using known values for the coefficients {a n ,b n ,n = 0,N — 1} 
for N = 1,2,3. Now knowing 2N coefficients in total allows us to estimate the next iV 
coefficients for either the a n or b n series. Table I shows the 'predicted' values versus the 
known values for N = 1. For the vector state with N — 1, the predicted values are in quite 
good agreement with the known values, to within a few percent. At iV = 2 (see Table II), 
however, a large discrepancy becomes evident between the predicted and 'known' values for 
the coefficient b<i - even the sign is different. This prompted us to re-examine the calculation 
of the weak-coupling coefficients b n by Sriganesh et al. [11]; and indeed we discovered some 
errors in those calculations. A corrected version of the calculations is given in the Appendix; 
the corrected series coefficients are also shown in Table II. 

For the scalar state, the 'predicted' series coefficients disagree with the known ones even 
at N — 1 - but this is easily understood, in that the simplest interpolation between the 
zeroth order series at either end would have no 'hump' in the middle such as is shown in Fig. 
4. The predicted values at iV = 2 agree qualitatively with the (corrected) known values. 

The Feynman-Kleinert approximants were found to give smooth and apparently accurate 
interpolations of the energy eigenvalues at finite m/g, up to order N = 2. At order N = 3, 
no solution was found to the system of simultaneous equations - we have not explored in 
detail why this occurred. Figures 3 and 4 show the results for the vector and scalar state 
interpolations. For the vector state energies, our interpolation is compared to data obtained 



-) 

mj 



2/3 



(24) 



6 



from the density matrix renormalization group (DMRG) method [8], the "fast moving frame" 
estimates of Kroger and Scheu [12], and the "renormal-ordered" mass perturbation series 
results of Adam [14]. The scalar state interpolations are compared to the finite lattice 
calculations of Sriganesh et al. [11], and again results of Kroger and Scheu, and Adam. Both 
plots show excellent agreement between the interpolation and the numerical data, within 
errors. We also plot the second order strong and weak coupling series used to generate the 
interpolation on the same plots, which diverge away for large and small masses respectively. 
The "renormal-ordered" expansion [14] about m/g = 0, extends the region of validity of the 
perturbation theory, and agrees fairly well with our interpolation for small masses. 

The interpolation performs particularly well for the scalar state, where for small m/g, 
the original strong coupling series disagrees significantly with the numerical data. The 
interpolation on the other hand agrees extremely well with existing results, despite the large 
gap that must be bridged between the two series. Tables III and IV list numerical values for 
specific values of m/g for the interpolation compared to previous works. We estimate the 
error on our interpolations based on the difference between our quoted N = 2 calculation 
and the corresponding N — 1 calculation. Error estimates have probably been overestimated 
somewhat using this method, as can be seen from the good agreement with the numerical 
data in Figure 4. We have not however found a better way of estimating these errors. 

IV. CONCLUSIONS 

The Feynman-Kleinert approximants have earned a mixed scorecard in the applications 
we have discussed. For the extrapolation of lattice strong-coupling series to the weak- 
coupling continuum limit, the approximants have basically failed. They converged for the 
ground-state eigenvalue, but only slowly; and for the first excited state eigenvalue, they did 
not converge at all. This is probably due to the presence of exponentially decaying terms 
corresponding to essential singularities in the weak-coupling limit, which are not accounted 
for in the Feynman-Kleinert approach. The upshot is, unfortunately, that this approach does 
not seem to have any advantage over previous techniques [7], which used Pade or integral 
approximants to extrapolate to intermediate couplings, and then 'matched' the results onto 
a weak coupling form to reach the continuum limit. 

For the continuum field theory, however, the Feynman-Kleinert approximants have 
proved useful in providing a sensible and apparently accurate interpolation between the 
known weak coupling and strong coupling series. The fact that we were able to pick up 
a mistake in the series coefficients by this approach adds confidence in the technique. No 
doubt many applications of a similar nature are possible. 
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APPENDIX A: WEAK COUPLING EXPANSION OF THE MASSIVE 

SCHWINGER MODEL 



Our starting point is the Hamiltonian as derived by Coleman in the two-particle subspace 
(see Eq. (4.18) of Ref. [10]). We set the background field 9 = 0. Expanding this Hamiltonian 
for large mass m, and keeping terms to order 0(m~ 3 ) we obtain a Schrodinger equation 



where 



{H + H 1 )^(x)=EiP(x), (Al) 



1 d 2 q 2 

F ' = -ma? + 7M (A2) 

mn Am 6 ox* Am z Zirm 6 ox z 



Rescale to dimensionless variables as follows: 



Then (Al) becomes 



X (mg 2 ) Z E (4m) 6 ' 



(H > + H[)^(z) = e^(z), (A5) 



where 



K = ~ + \z\ (A6) 



The last term in (A2) has now been dropped as it is now clear it is of order ~ (g/m) 2 . The 
solution of the leading-order Schrodinger equation 

" m) 1>(z) = -eiP(z), (A8) 

was discussed by Hamer [13]. The solution of (A8) is a symmetric or antisymmetric Airy 
function which obeys the condition 

Ai'(-e„) = (symmetric) (A9) 
Ai(— e n ) = (antisymmetric). (A10) 

The lowest symmetric (antisymmetric) wavefunction gives the energy for the vector (scalar) 
state. 

The higher order terms contained in H[ may be taken into account by calculating the 
correction term Ae = (ifj\H[\ilj) using numerical integration (the results are shown in Table 
V). Rescaling back to our original variables gives us 
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= 0.6418 + . 15 47 (All) 



for the vector state binding energy, and 

1 (9 \ nmno f 9 

V 



* = 1.473 (i-f-lfJi) _ 0.1093 m 5/3 . (A12) 
g \mj 7T \mj \mj 



for the scalar state binding energy. 
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TABLES 



TABLE I. Comparison of known coefficients in weak [11] and strong [9] coupling series, versus 
the Feynman-Kleinert estimates (labeled as "Predicted"). Feynman-Kleinert estimates are ob- 
tained by starting with zeroth order coefficients ao and bo (i.e. N = 1), for each of the vector and 
scalar particle states. 



Coefficient 




Vector State 


Scalar State 




Predicted 


Known 


Predicted 


Known 


ai 


-0.2046 


-0.2189 


-0.2709 


1.562 


bi 


-0.2917 


-0.3183 


-0.8813 


-0.3183 


TABLE II. 


As for Table I 


, but with N = 2. The "known" coefficients for 62 


contain both the 


result obtained 


in Ref. [11] (in square brackets), and our 


recalculation (unbracketed) . 


Coefficient 




Vector State 


Scalar State 




Predicted 




Known 


Predicted 


Known 




0.2298 




0.1907 


-10.13 


-13.51 


03 


-0.3148 






49.01 




b 2 


0.2159 




0.1547 [-0.2521] 


-0.1854 


-0.1093 [0.1085] 


b 3 


-0.1669 






0.2794 




TABLE III. 


Comparison of interpolation results for the vector state with numerical results 


obtained by Byrnes et al. [8], 


Kroger 


and Scheu [12], and Adam [14]. 




m/g 


This work 




Byrnes 


Kroger and 


Adam [14] 








et al. [8] 


Scheu [12] 




0.125 


0.540(2) 




0.53950(7) 


0.528 


0.539 


0.25 


0.520(3) 




0.51918(5) 


0.511 


0.516 


0.5 


0.490(4) 




0.48747(2) 


0.489 


0.474 


1 


0.448(4) 




0.4444(1) 


0.455 


0.43 


2 


0.396(4) 




0.398(1) 


0.394 


0.49 


4 


0.341(3) 




0.340(1) 


0.339 


0.76 


8 


0.287(2) 




0.287(8) 


0.285 


1.40 


16 


0.237(2) 




0.238(5) 


0.235 


2.75 
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TABLE IV. Comparison of interpolation results for the scalar state with numerical results 
obtained by Sriganesh et al. [11], Kroger and Scheu [12], and Adam [14]. 



m/g 


This work 


Sriganesh 
et al. [11] 


Kroger and 
Scheu [12] 


Adam [14] 


0.125 


1.23(13) 


1.22(2) 


1.314 


1.220 


0.25 


1.24(17) 


1.24(3) 


1.279 


1.230 


0.5 


1.21(19) 


1.20(3) 


1.227 


1.165 


1 


1.12(17) 


1.12(3) 


1.128 


0.99 


2 


0.99(13) 


1.00(2) 


0.991 


0.88 


4 


0.84(8) 


0.85(2) 


0.837 


1.07 


8 


0.69(5) 


0.68(1) 


0.690 


1.78 


16 


0.56(3) 


0.56(1) 


0.559 


3.41 



TABLE V. Results of numerical integration over symmetric and antisymmetric states for two 
operators. 



Integral 


Symmetric 


Antisymmetric 




-0.577655 


1.093349 


<*(*)>o 


0.490777 


0.0 



12 



FIGURES 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 

1/X 

FIG. 1. Ground state energy per site of the lattice Schwinger model at m/g = 0, as obtained by 
the Kleinert approximants with N = 5 and N = 13. Open circles denote finite lattice calculations, 
the closed circle shows the position of the exact result in the continuum limit, uq/xN = —2/ir. We 
also show an extrapolation using a standard Pade approximant for comparison. 
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/ 



-X 



N=3 
N=4 
N=5 
N=6 
Exact 



0.9 



1.1 



1.2 



1.3 



1.4 



1.5 



1.6 



1.7 



1.8 



FIG. 2. Dependence of the value &o as a function of Co for the case m/g = 0, at four orders in 
N. The exact result is shown as a solid line at bo = ~ 1.128. 
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t 1 1 ; r 



/g 




0.2 " 



-3-2-101234 

log 2 (m/g) 

FIG. 3. A comparison of the Feynman-Kleinert interpolation for the vector state with numerical 
DMRG data of Byrnes et al. [8], "fast moving frame" estimates of Kroger and Scheu [12], and 
the "renormal-ordered" mass perturbation theory results of Adam [14]. The 2nd order (N = 2) 
weak and strong coupling series, which are used to generate the interpolation, are also shown for 
comparison. 
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1.6 



1.4 - 




/ 2nd order weak coupling expansion 



-3-2-10 1 2 3 

log 2 (m/g) 

FIG. 4. A comparison of the Feynman-Kleinert interpolation for the scalar state to finite lattice 
estimates of Sriganesh et al. [11], "fast moving frame" estimates of Kroger and Scheu [12], and 
the "renormal-ordered" mass perturbation theory results of Adam [14]. The 2nd order {N = 2) 
weak and strong coupling series, which are used to generate the interpolation, are also shown for 
comparison. 
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